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Compact binaries consisting of neutron stars / black holes on eccentric orbit undergo a perturbed 
Keplerian motion. The perturbations are either of relativistic origin or are related to the spin, mass 
quadrupole and magnetic dipole moments of the binary components. The post-Newtonian motion of 
such systems decouples into radial and angular parts. We present here for the first time the radial 
motion of such a binary encoded in a generalized Kepler equation, with the inclusion of all above- 
mentioned contributions, up to linear order in the perturbations. Together with suitably introduced 
parametrizations, the radial motion is solved completely. 



I. INTRODUCTION 

The worldwide effort to capture gravitational wave sig- 
nals emitted by astrophysical sources is under way. A 
network of interferometric gravitational wave detectors 
0-0 is either already operational, or close to comple- 
tion. Compact binaries consisting on neutron stars / 
black holes are among the best candidates to emit grav- 
itational radiation in the bandwidth of these detectors. 
Upper limits on the gravitational radiation emitted by 
such binaries were already found , from the S2 sci- 
entific run of LIGO. 

The evolution of such a binary system can be divided 
into three phases: inspiral, merger and ringdown. The 
merger phase can be understood only by numeric simu- 
lations. Even in the last part of the inspiral a numerical 
treatment seems adequate due to the intermediate bi- 
nary black hole (IBBH) problem 0, studied also in 81. 
Earlier, in the inspiral phase a post-Newtonian (PN) de- 
scription of high accuracy provides satisfactory results. 
The PN corrections of relativistic nature are known to 
3PN orders . However there are other contributions to 
be taken into account, related to various physical char- 
acteristics of the binary components as well. 

For compact binaries, there is the spin-orbit (SO) in- 
teraction appearing at 1.5 PN orders. At this accuracy 
both the spin vectors Si and the orbital angular momen- 
tum L undergo a precessional motion about the total 
angular momentum J [lfj. This is a novel feature in the 
post-Newtonian evolution of the system. Such an effect 
was recently claimed [ill to be observable for the J0737- 
3039 A/B double pulsarQ, 0. 

The precessional motion of the spin(s) is called simple 
precession, whenever the two masses are equal, or one of 
the spins can be neglected, say S2 = These two 
cases were studied in [151 . where among other results, 
the Kepler equation was derived up to 3PN orders with 
the inclusion of the SO contributions. As a related re- 
sult, the evolution of the relativistic periastron advance 
parameter was recently computed [Ig. The tilt angle of 
the spin with respect to L was estimated to be smaller 
than ~ 60° from generic astrophysical considerations on 
the evolution of compact binaries |l7j . 

The SO interaction gives corrections to the losses of en- 



ergy and magnitude of angular momentum of the system 
occurring due to gravitational radiation. For eccentric 
orbits these were given by an d An other work 
relying on the use of the Effective One-Body approach 
[2Cj . |2l|. has employed the SO contribution in the study 
of the inspiral to plunge phase of the coalescence . 

Moreover, in the two cases of simple precession de- 
tection template families have been worked out both 
containing a set of phenomenological parameters |2^| or 
physical parameters j24[. The latter would allow for de- 
termining the angle k,\ and the magnitude of the single 
spin (in fact of x — Si/m\) from the study of gravita- 
tional radiation. 

At 2PN an other set of new effects related to various 
physical characteristics of the compact binary emerge. 
The losses of energy and magnitude of angular momen- 
tum of the system on eccentric orbit, due to gravitational 
radiation were derived in |25| for the spin-spin interac- 
tion, in [2^| for the mass quadrupole-monopole interac- 
tion and in [2]j for the interaction of magnetic dipoles. 

Thus at 2PN physical quantities like the mass 
quadrupole and magnetic dipole moments, as well as an- 
gular variables characterizing the spins and moments ap- 
pear in the formalism. Neither detection templates, nor 
methods to find out these new physical parameters have 
been worked out so far. We note that in principle, the 
observation of the evolution of gravitational wave fre- 
quency |28| allows to impose constraints on a combina- 
tion of these parameters, but does not allow to predict 
their individual values. There is still much to do until 
a complete understanding of the complicated motion the 
system, occurring when all these interactions are taken 
into account, will be achieved. Our present work fills 
an important gap in the description of compact binaries 
with the enlisted physical characteristics. 

In Section 2 we describe the post-Newtonian motion 
of such a binary system. The radial part of the motion 
decouples and defines a radial orbit. We give the generic 
expressions of the turning points of the radial motion 
which allows for the introduction of the semimajor axis 
a r and radial eccentricity e r . We define the generalized 
true anomaly \ arid eccentric anomaly £ parameters. Our 
generalized true anomaly parameter is different from the 
one employed in the Damour-Deruelle formalism and we 
establish their relation. In Section 3 we derive the main 
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result of the paper, which is the generalized Kepler equa- 
tion 



n(t-t ) = f — e t sin£ 

+/ t sin + -?)] 

2 

+ Y, ft sin [ X + 2(^o-^)] • (1) 

i=l 

The explicit expressions of the orbital elements n, et, ft 
and f\ are given in Section 4, together with all contribu- 
tions to a r and e r . They contain the relativistic contri- 
butions (PN-terms), together with the SO, SS, DD and 
QM terms, all to linear order. One of the consequences of 
taking into account the physical characteristics of the bi- 
nary, like the spins, mass quadrupole and magnetic dipole 
moments is the emergence of novel angle variables in the 
Kepler equation QJ. The angles ipi are the azimuthal 
angles of the spins and 2tp = ipi + ifa. (For more details 
on the notations see The angle ipo is the argument 

of the periastron (defined here as the angle subtended 
by the periastron and the intersection line of the planes 
perpendicular to the total and orbital angular momenta, 
respectively) . 

We emphasize that in the description of the SO interac- 
tion several spin-supplementary conditions (SSC) can be 
used. In order to simplify the formalism, in this paper 
we use the non-covariant SSC of Pryce |29l ] and New- 
ton and Wigner [3(j ■ In this SSC the Lagrangian is not 
acceleration-dependent and the radial equation is simpler 
than in the covariant SSC p9| employed earlier in 



II. GENERALIZED TRUE AND ECCENTRIC 
ANOMALY PARAMETRIZATIONS 

The linear contributions to the motion of the compact 
binary can be collected in the Lagrangian 
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Note that the SO-part of the Lagrangian above was not 
given before and it is valid when the spin-supplementary 
condition (SSC) of Pryce and Newton and Wigner 
|30| is chosen. We have verified that the SO-part of the 
acceleration derived from Cso agrees with the expres- 
sion (Alb) of [23|. The magnitude and direction of the 
spins are denoted as Si and Si. The angle subtended 
by them is 7 = cos _1 (Si ■ S 2 ). Here S = Si + S 2 
and a = (777,2/777.1) Si + (7711/7772) S 2 . The magnitude 
and direction of the magnetic dipole moments di are 
denoted as di and di. They subtend the angle A = 
cos _1 (di • d 2 ) with each other. In a coordinate systems 
K, with the axes (c, L x c, L), where c is the unit vector 
in the J x L direction, the polar angles Ki and tpi of the 
spins are defined as Si =(sin Ki cos ipi, sin Ki sin ip^, cos Ki) 
(see In tn e coordinate system K % with the 

axes (bi, Si x bi, Si), where bi arc the unit vectors 
in the Si x L directions, respectively, the polar angles 
cti and Pi of the the magnetic dipole moments di are 
d ; =(sinajCos/3j,sinQ!jSin/3j,cosaj) (see [13 )■ The 
quadrupolar parameters (see |26[l are defined as Pi = 
Qi/rriim 2 , where Qi is the quadrupole-moment scalar [33] 
of the i th axially symmetric binary component with sym- 
metry axis Si. The reduced mass is [i = mxmilvn and 
77 = /i/777. 

From a radial equation can be derived 
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Here L — (l/2w) L(x)dx is the angular average of the 
magnitude of orbital angular momentum x being 

the true anomaly parameter. The explicit values of L in 
the case of spin-spin, quadrupole-monopole and magnetic 
dipole-dipole interactions were computed in |2J| , [26] and 
[23 • A is the magnitude of the Laplace-Runge-Lenz vec- 
tor characterizing a Keplerian motion with E and L. The 
coefficients SAi in Eq. are constant PN perturbations 
given in the Table which can be read from [3lT | . The 



TABLE I: Various post-Newtonian constants in SAi. 



5A 


3(3r7-l)^ 


5A 1 


2(777 - 6)^ 


SA 2 


-2(3r? + (577 - 10)^1 


SA 3 


(-3r7 + 8)^^ 



SO, SS, QM and DD contributions to SL and SE are en- 
listed in the Tables ITT1 and HTH The shorthand notations 
oidd and /3dd(2x) are defined in Eq. 10. 

In |3(j a generic scheme was introduced for parametriz- 
ing such perturbed Keplerian motions. The advantage of 
the generalized eccentric anomaly parametrization and 
generalized true anomaly parametrization is that a sim- 
ple technique based on the residue theorem can be ap- 
plied for computing secular effects j3|| . The scheme was 
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TABLE II: Various contributions to SL. 



SO 


— 2 

G^L l^-^-Si COS Ki 


ss 


-^jS^iSinm sinK 2 {2Acos [ X + 2 (i>o 
+ (3Gm.fi + 2A cos x) cos 2 ( X + 0o - 4>) } 


QM 


ELiP^in 2 Ki {2Acos [x + 2 (jfa - ^)] 
+ (SGmfj, + 2y4cosx) cos 2 (x + Vto — 4>i)} 


DD 


fi 2 d 1 d 2 
2L a 


(3Gm^+4A cos x )Pdd (2 X ) - A sin X dfe jf x) 



TABLE III: Various contributions to 8E. 



SO 


no contribution 


ss 


— ^i^i { 3 COS Kl COS K2 — cos 7 

—3 sin fti sin K2 cos 2 (% + 0o — W) } 


QM 


^r-ELiW [1 -3sin 2 Kl cos 2 fr+ 


DD 


4j# [a OD - 3/3 DD (2*)] 



applied individually for each of the SO, SS, QM, DD per- 
turbations in 0, [25j, [2(| and [23, respectively. It is 
straightforward to derive the PN contribution to these 
parametrizations. We give here in a concise form both 
parametrizations, with the inclusion of the PN contri- 
bution as well. Both are defined in terms of the turning 
points r m j n and r max of the radial motion, given by r = 0: 
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where 

Pi = sin at cos Pi , 

(Tj = cos a, sin Kj + sin ctj sin /3j cos . (7) 

The generalized eccentric anomaly parametrization r (£) 
is then defined as 



r (£) = a r (1 — e r cos£) 



(8) 



The eccentric anomaly £ reduces to the eccentric anomaly 
parameter u of j^l for 1PN perturbations. In Eq. © 
the semimajor axis a r and the radial eccentricity e r was 
introduced as 
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The generalized true anomaly parametrization is defined 
as: 



r(x) 
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These two parametrizations of the radial motion are in- 
terrelated by the Keplerian relations 
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with e r in place of the Keplerian eccentricity. Note that 
in the Damour-Deruelle formalism [3jJ , a different gener- 
alized true anomaly parameter v is introduced. To 1PN 
accuracy v is related to x as 



, X 
tan — = 
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' Gmn / 2 3 12EL^ 
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When using the generalized true anomaly parameter v, 
the equation £ (v) replacing Eq. I|12|l will contain the 
angular eccentricity ee rather than e r . 



We have introduced the shorthand notations 

a SS = 3 COS K\ cos Ki — cos 7 j 

flss = sin^! sinK2 cos 2 (ip — , 
a QM = 2-3 sin 2 Ki , 
@q M = sin 2 m cos 2 (-00 - tpi) , 
ocdd = 2 cos A + 3(pi(72 — P2C1) sin 
— 3(pi/?2 + cos Atjj , 
PDD(kx) = {o\(J2 - pxp-z) cos [k\ + 2 (ipo - V 7 )] 
-(pi&2 + P2<Ji) sin [fcx + 2 (^o - V>)] 
= /3co(0) , 



III. GENERALIZED KEPLER EQUATION 

The first three terms on the right hand side in the 
radial equation (@J sum up to: 



2E 2Gm 



L 2 



2 ™2 



2.4 



1 r a* / m 



Asin 2 x (<5Q+^Pcosx) ,(15) 



' where <$P and <5Q are perturbation terms depending on 
(6) the physical parameters of the binary and they are given 
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TABLE IV: Various contributions to SP. 



PN 



SO 



ss 



QM 



DD 



A (r) - 2) - 4G 2 mV 
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_ y 



Si cos Ki [A +3G 2 m 2 fj,' 



G^ i S 1 S 2 
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+/%mU 2 + G 2 ™V 
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TABLE V: Various contributions to 



PN 


--^[2G 2 mV^(3r,-19) 
-G 4 mV (77 + 9)+ A 4 (37? -1)1 


SO 


2^f%. S m , 3 5,cos^( 

8=1, J^I 


'3A 2 + G 2 mV) 


SS 


G 2 m,t 4 SiS 2 
2c 2 r B A 


a ss (3A 2 +G 2 mV 


) + 2f3 QM A 2 




QM 


G 2 m 4 ,i 5 V^ 2 


ajjM (3A 2 + G 2 mV) + 2/3 ss A 2 


DD 


2L 6 A 




;3A 2 + G 2 mV 


2 )+2f] DD A 2 





in the Tables llVl and Ivl The parametrization (jl 1|> has 
the advantageous property 



dr 
dx 



r smx 



(16) 



Employing Eqs. (|11() and (|16|l into the radial equation 
Q , then taking the square root and forming the recipro- 
cal, after a series expansion to first order in the pertur- 
bations we find 



dt 

dx 



L \ 2\i 2 r 2 A sin 2 x 

3 

-2 (L5L + r 2 ^iSE) + ^ 5 A 



r 2 nA(8Q+6P cosx) 



(17) 



The integration of this equation is the main purpose of 
our paper. 

Apparently Eq. (|17fl becomes singular at X — k7T, 
k G Z , because of the sin 2 x term in the denominator. 
Such singularities could be just apparent, as was shown 
for the SO-contribution in 36], We have verified that 
after forming the common denominator in the bracket, 
and inserting the detailed expressions of 6P, SQ, SL and 
SE, given in Tables IIIIVI the numerator becomes pro- 
portional to sin 2 x- Therefore the singularities are just 
apparent rather than real for all type of contributions 
considered here. Thus we obtain 



dt _ fir 2 ( dt\ / ' dt\ ( dt\ ( dt \ / dt \ 

dx L \dx)p N \dx)so \ d x) S s \ d x) DD \dx) QM 

dt\ = j»_ , -13)G 2 mV + (3 r?-l)A 2 + ( 3r? - 8) Gm/M cos x I , 
dxJ pn 2c 2 L 



( dt\ G[i 2 r 2 . -r s ^ Ami + 3m., 

— = o- (3Gmu, + A cos x) > -biCOSKi , 

G/x 3 5i5 2 r 2 f r ^ .„ ^ . — j- . „ . n 2l 4 . . 1 

I — J = =g — < |_Gm/i(3ass + Wss) + A (a S s + Pss) cos x\ ^TJ sm/ci sihk 2 cos|_x + 2 (V'o ~ i>)\ > , 

\ a X/ ss 2c 2 L [ ii r A J 

f dt\ Grn^i^r 2 2 f — 2L 4 ] 

— J = =5 — 22Pi{[ Gm K 3a QM+^QM) +A(a t QM +0 t Q lu )coBx] - -f^=F sin 2 cos[ X + 2 (V'o -*l>i)] > , 

\ a XJ qm 4L i=1 { fi r A J 

/ dt\ = _ /i did 2 r | (3a ^^ + + ^ + ^ + cos [ X + 2 (Vo - ?)] ) ■ (18) 



These expressions are regular for any x- the paper. 

By integrating Eq. (|18|l and employing the relation 
(|12|l between the two parametrizations we obtain the gen- 
eralized Kepler equation ft). This is the main result of 
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IV. THE ORBITAL PARAMETERS 



In this Section we enlist the detailed expressions of 
the orbital parameters appearing either in the eccentric 
anomaly parametrization J5J) or in the generalized Kepler 
equation 0). 

The semimajor axis is 
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The radial eccentricity is 
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(22) 



In what follows, we enlist the parameters / t and f\, 
the analogues of which also appear in the extension to 
2PN j3?J of the Damour-Deruelle parametrization, how- 
ever have no PN counterpart. The parameter f t receives 
only SS and DD type contributions 



ft = ft S + f, 



cDD 



rSS 
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-2E\ 3/2 /ISA 



c 2 mAL 
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(23) 



Finally the parameters fl, originating from the QM in- 
teraction, are given as 



ft 



2S\ 3/2 toV 



PiSin Ki 



(24) 



CONCLUDING REMARKS 
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The mean motion is 
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1 - (?7 - 15) 
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The generalized Kepler equation with the or- 
bital parameters (|19fl - H24fl . together with any of the 
parametrizations JSJ or l|ll(l and their relation (|12(l rep- 
resent the complete solution of the radial motion of the 
compact binary on eccentric orbit, to linear order in the 
perturbations. All perturbations arising from relativis- 
tic corrections and from the presence of spins, mass 
quadrupole and magnetic dipole moments are included 
here to linear order (PN, SO, SS, QM, DD contributions). 

The generalized Kepler equation contains two parame- 
ters, the generalized eccentric anomaly £ (this is defined 
similarly as the parameter u in |3lL 1371 ] and the general- 
ized true anomaly \ (different from the parameter v of 
[3TI Wfy. their relation being given to 1PN accuracy by 
Eq. (Hi). 

The generalization of the relation l|14|) to 2PN accu- 
racy would imply to give the \~ parametrization Ijlll) to 
2PN, however by the method described in |3jj r (x) can 
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be defined only to linear order in the perturbations. Nev- 
ertheless, the linear contributions to the Kepler equation 
included in the present paper (containing the parame- 
ters £ and x)i an d the 2PN contributions [33 (containing 
u = £ and v) can be simply summed up, at the price of 
having all three parameters (£ = u, x an d v) present in 
the formalism. Notably at 2PN orders a new parameter, 
gt is also present (3^1 • We note here that the 3PN contri- 
bution to the Kepler equation is also known [l5j , however 
the contribution of the first PN correction of the SO in- 
teraction, arising at 2.5 PN is not. Therefore we conclude 
that for the moment, the radial motion is solved only to 
2PN orders accuracy. 

In contrast with the PN and 2PN Kepler equations, our 
Eq. contains the additional angles ipo and ipi. Du ring 
one radial period these can be considered constants [19j. 
On the long run however, all these angles slowly vary 
as the orbit and the spins undergo precessional motions. 
In order to describe the slow evolution of these angles, 
the study of the angular part of the motion (as opposed 
to the radial one) is necessary, with the inclusion of all 
perturbations, to linear order. This is available for the 
PN perturbation [3lJ and SO perturbation [l^, however 



the latter holds only for special cases (equal masses or 
a single spin). A systematic investigation of the angular 
part of the perturbed Keplerian motion is under way • 
We remark that the circular orbit limit of our formu- 
lae should not be taken as A — ► 0. This is because the 
standard interpretation of the Laplace-Runge-Lenz vec- 
tor holds only in the Newtonian limit. In the circular 
orbit limit all corrections considered here add nonvanish- 
ing contributions to A. 

Our generic Kepler equation and the orbital el- 
ements (I19ll - H24|) correctly reproduce the PN contribu- 
tions |3lj| and the SO contributions 0]. However all of 
the other contributions are new. 
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